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Abstract 

Hill's equations arise in a wide variety of physical problems, and are specified by 
a natural frequency, a periodic forcing function, and a forcing strength parameter. 
This classic problem is generalized here in two ways: [A] to Random Hill's equations 
which allow the forcing strength q^, the oscillation frequency and the period 
(At)^ of the forcing function to vary from cycle to cycle, and [B] to Stochastic Hill's 
equations which contain (at least) one additional term that is a stochastic process 
£. This paper considers both random and stochastic Hill's equations with small 
parameter variations, so that Pk = Qk~ (<Zfc)> = (Afc), and £ are all 0(e), where 
e<l. We show that random Hill's equations and stochastic Hill's equations have the 
same growth rates when the parameter variations pk and Ik obey certain constraints 
given in terms of the moments of £. For random Hill's equations, the growth rates 
for the solutions are given by the growth rates of a matrix transformation, under 
matrix multiplication, where the matrix elements vary from cycle to cycle. Unlike 
classic Hill's equations where the parameter space (the X-q plane) displays bands of 
stable solutions interlaced with bands of unstable solutions, random Hill's equations 
are generically unstable. We find analytic approximations for the growth rates of 
the instability; for the regime where Hill's equation is classically stable, and the 
parameter variations are small, the growth rate 7 = 0(p 2 ) = 0(e 2 ). Using the 
relationship between the [&kiPk) an d the £, this result for 7 can be used to find 
growth rates for stochastic Hill's equations. 
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I. INTRODUCTION 

Hill's equation is a second order periodic differential equation that arises in many 
physical problems [22]. In addition to its original application for lunar orbits [15], Hill's 
equation describes celestial dynamics [7,8,21], orbit instabilities in dark matter halos [5], 
particle production at the end of the inflationary epoch in the early universe [17,18,19], 
the motion of a jogger's ponytail [16], and many other physical systems. The goal of this 
paper is to consider the addition of random perturbations to Hill's equation and thereby 
generalize our current understanding. Toward this end, we adopt the following definitions: 

Definition 1.1: A Classical Hill's equation has the form 

+ [A + qQ(t)]y = , (1) 

where the parameters (A, q) are constant and where the barrier shape function Q{t) is 
periodic, so that Q(t + At) = Q(t), where At is the period. For the sake of definiteness 
we take At = n and normalize the barrier shape function so that Jq Q{t)dt = 1. 

Definition 1.2: A Hill's Equation is said to be Cycle-to- Cycle Random (CC Random) 
if the differential equation is deterministic over each interval of the periodic function Q, 
but the values of the forcing strength q k and/or oscillation frequency X k take on different 
values for each periodic interval. The CC Random Hill's equation has the form 

d ^ + [\k + qkQ(t)]y = 0, (2) 

where the parameters (Afc, q k ) vary from cycle to cycle, and are drawn independently from 
well-defined distributions. The index k labels the cycle. 

In this work the period At is considered fixed; one can show that cycle to cycle varia- 
tions in At can be scaled out of the problem and included in the variations of the (X k , q k ) 
by changing their distributions accordingly (see Theorem 1 of Ref. [2]). 

Periodic differential equations in this class can be described by a discrete mapping of 
the coefficients of the principal solutions from one cycle to the next. The transformation 
matrix takes the form 

h k {h 2 k - l )/9k 
.9k h k 

where the subscript denotes the cycle. The matrix elements for the kth cycle are given by 



(3) 



h k = 2/i(tt) and g k = yi(n) , (4) 

where y\ and y2 are the principal solutions for that cycle. The index k indicates that the 
quantities (\k,qk), and hence the solutions (h k ,gk), vary from cycle to cycle. Throughout 
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this work, the random variables are taken to be independent and identically distributed 
(iid). 

Note that the matrix in equation ([3]) has only two independent elements (not four). 
Since the Wronskian of the original differential equation (J2J) is unity, the determinant of 
the matrix map ([3]) must also be unity, and this constraint eliminates one independent 
element. In addition, this paper specializes to the case where the periodic functions Q(t) 
are symmetric about the midpoint of the period. This property implies that 2/1(71") = y<i (ir) 
= hk, which eliminates a second independent matrix element [1,22]. These two constraints 
imply that 2/2 C 71 ") = (hi — l)/gk, resulting in the form for the matrix given by equation fl3]). 

The growth rates for Hill's equation (|2]) are determined by the growth rates for matrix 
multiplication of the matrices M. k given by equation (|3]). Here we denote the product of 
N such matrices as Ai^ N \ and the growth rate 7 is defined by 

7 = lim Ilog||MW||. (5) 

This result is independent of the choice of the norm 1 1 • 1 1 , and the limit exists almost surely, 
as shown in previous work [12,13,20]. 

Definition 1.3: A Stochastic Hill's Equation is a Hill's equation with the inclusion of one 
or more additional terms that is a stochastic process. In this paper, the stochastic Hill's 
equation is written as a Langevin equation. For example, one can write stochastic Hill's 
equations in the form 

+[A + gQ(t)h/ = £, (6) 

or 

^ + [A+(g + e)Q(t)]y = ) (7) 

where £ is a stochastic process and (A, q) are constant. 

Note that in the stochastic Hill's equation, we can evaluate the stochastic term £ using 
the calculus of either Stratonovich or Ito. For purposes of this work, however, we only 
require that the process £ has zero mean and finite variance. 

The goal of this paper is to study both the CC Random Hill's equation and the Stochas- 
tic Hill's equation in the limit where the parameter variations and/or the stochastic terms 
are small. For the CC Random Hill's equation, the parameters are rewritten in the form 

A fc = A + 4 and q k = q + p k , (8) 

so that (A, q) are constant whereas the l k and p k are small and are allowed to vary from 
cycle to cycle. The distributions of the i k and the p k are assumed to have zero mean 
= = (p k ) and finite second moments (£%) and {pi). 
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For both the CC Random and the stochastic Hill's equation, the solution over a given 
cycle (period of the function Q(t)) can be expanded such that 

y(t) = y {t)+y k (t), (9) 

where yo obeys the Classical Hill's equation (pQ), and is the same for every cycle. The 
correction term yk{t) is, in general, different for every cycle, due to the parameter variations 
(£k,Pk) for the CC Random Hill's equation and due to different realizations of the stochastic 
process £ for the Stochastic Hill's equation (different samplings of the process). 
Next we introduce an ordering parameter e C 1 such that 

y k = 0(e), 4 = 0(e), Pk = 0(e), and £ = 0(e) . (10) 

The results of this paper are correct to leading order in e. 

Remark: The subscripts '1' and '2' on the function y(t) denote the two principal solutions 
of any Hill's equation for a given cycle, i.e., for given values of the parameters. The 
subscripts '0' and on y(t) denote the leading order and first order parts (in e) of the 
function, respectively (see equation (J9])). 

To place this work in context, we note that a number of previous studies have con- 
sidered Mathieu's equation, a particular type of Hill's equation, with stochastic forcing 
[25,26]. These papers show that stochastic fluctuations can lead to instability in otherwise 
stable systems. Other related work includes random operators in the the Schrodinger wave 
equation in quantum mechanics and population variations in biomathematics [6,10,23], and 
stability studies of the quasi-periodic Mathieu equation [28]. This body of existing work 
does not, in general, provide closed-form expressions for the growth rates of the differential 
equations. Our previous work has studied the CC Random Hill's equation in the regime 
where the forcing strengths are large and the solutions are highly unstable [2] , for the limit 
where the periodic functions Q(t) become (periodic) delta functions [3], and for more gen- 
eral cases [4]; here we find closed- form expressions for the growth rates of corresponding 
random matrices (given by equation (j3J)). Examples of such analytic results are rare in the 
literature - the classic papers [9,12,13] show the existence of growth rates, but present few 
examples (see also Ref. [24]). 

The main goal of this present work is to build a bridge between the Stochastic Hill's 
equation and the CC Random Hill's equation. This paper focuses on the regime of small 
fluctuations as defined through equation (flOl) . We derive the conditions necessary for 
the two types of Hill's equations to have the same growth rates (Theorem 2.3). We also 
find analytic expressions for the matrix elements (Theorem 2.1) and the growth rates 
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(Theorem 2.2) for CC Random Hill's equations. The combination of these three results 
provides analytic expressions for the growth rates for Stochastic Hill's equations. Section III 
presents a special case of the problem where the natural oscillation frequency A is constant 
and the forcing parameters q k are small but fluctuating from cycle to cycle. The paper 
concludes in Section IV with a summary and short discussion of the results; representative 
applications are described in the appendices. 

II. HILL'S EQUATION WITH SMALL FLUCTUATIONS 

In this section, we consider the case where the parameters of Hill's equation are arbi- 
trary, but the variations from cycle to cycle and/or the amplitude of the stochastic process 
are small. 

Theorem 2.1: Consider the CC Random Hill's equation (T5]) with the ordering of equation 
( jTUl) . Let h and g denote the matrix elements of the map (J3j), and let ho and go denote the 
matrix elements in the absence of fluctuations [l k = = pt). Provided that the zeroeth 
order matrix elements h and go are nonvanishing, the matrix elements h and g are given 
by 

h = h - i k X - p k Y + 0(e 2 ) and g = g - W - PkZ + 0(e 2 ) , (11) 
where X, Y, W, Z are integral quantities that are constant from cycle to cycle: 



X = , 

2g Jo 



W 



2 h Jo 



u2 i \ 2 2 2 

K~ l ) ym-%Vo2 
(K + l) Vol - 9oVo2 



dt . 



dt . 



1 f n - 
2g Jo 



U2 i \ 2 2 2 

K ~ 1 ) 2/oi - 9oVo2 



2ho Jo 



Q(t) [(hl + l) yl x -glyl 2 



dt , 

dt. 
(12) 



Proof: By definition, the leading order part of the solution yo obeys the Classical Hill's 
equation in the form of equation ([T]), where A and q are constant. To all orders (in e), and 
for both principle solutions, the correction y^ obeys a differential equation of the form 

d 2 y k 



dt 2 



+ 



X + £k + (q + Pk)Q 



Vk + 



+ PkQ yo = o. 



(13) 



where yo is the zeroeth order part of the principal solution. To leading order in e, this 
equation reduces to the form 

d 2 y k 



dt 2 



+ 



X + qQ 



yk + 



+ p k Q y + 0(e 2 ) = . 



(14) 



For the remainder of this argument, we neglect terms that are 0(e 2 ), so that the results 
are correct (only) to first order in e. After multiplying equation (|14p by yo and integrating 
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over one cycle, the expression becomes 
d 2 y k 



dt 2 



y dt + X + qQ yoykdt + 



+ PkQ vldt = 



(15) 



where we ignore higher order terms. By integrating the first term by parts, applying the 
boundary conditions at t = 0, and using the fact that yo satisfies the zeroeth order equation, 
we obtain 

yoMz/fcM - yo(n)yk(n) + 4 / vldt + p k Qvtdt = o . (16) 

Jo Jo 

This equation holds for both cases, i.e., where y = yo + yk is the first or the second principal 
solution. Using the definitions of the matrix elements, we obtain two equations for two 
unknowns (note that we suppress the index k as convenient to simplify the notation) 



and 



h 2 
"o 



ho9 - 9oh + K x = 



1, , h 2 -l T ^ 
—h - h h K 2 = 



where the integrals K\ and K 2 are defined by 



ylidt and K 2 



v -k+VkQ yl 2 dt 



The matrix element h is thus given by the quadratic equation 



h 2 + h 



9o 



K 1 - g Q K 2 



+ K X K 2 - h 



o = 0. 



(17) 
(18) 

(19) 
(20) 



which has solution 



h = ± 



hi 



4\ 9o 



-K x - g K 2 



1/2 



1 

2 



h ^K x - g K 2 
. 9o 



(21) 



We choose the + solution so that h — > h in the limit £k,Pk 0. This expression can be 
simplified by keeping only the leading order terms in e, i.e., 

1 



A similar result can be obtained for g, i.e. 



-K x - g K 2 + 0(e 2 ) . 



9 = 9o 



1 

2hv, 



{h 2 + l^-g 2 ^ +0{e 2 ). 



(22) 



(23) 



Since we want to isolate the dependence of the results on the parameters tk and pk that 
vary from cycle to cycle, the matrix elements can be rewritten to take the forms given in 
equations ffTTl) and f|T2l . • 
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Keep in mind that only the (£k,Pk) vary from cycle to cycle. As a result, the parameters 
X, Y, W, Z are constant, i.e., the same for all cycles. 
Alternatively, we can define the moment integrals 



yljdt and Jj 



Q(t)y 2 0j dt } 



(24) 



where j = 1,2, and can write the parameters X, Y, W, Z in the form 



X = 

2g 



hi 



1/1-^2 



Y 



1 

2^ 



(/ig-l) h-glh 



W 



1 

2h 



(h 2 + l)h- glh 



and 



1 

2h~o 



hr 



l)Jx-gp2 



(25) 



Remark: Note that this perturbation scheme fails in the limit where the leading order 
principal solutions vanish, i.e., when either ho — >* or g — > 0. 

The following results consider the case where \ho\ < 1 so that the solutions are classically 
stable (i.e., stable in the absence of fluctuations in the parameter values). For the case of 
nonzero parameter fluctuations, however, the solutions grow exponentially. 

Theorem 2.2: For a CC Random Hill's equation in the regime where the leading order 
matrix element \h \ < 1, the solutions are stable in the limit £ k ,p k — > 0. In general, the 
solutions are unstable with growth rate 



7 



ho 



1 2 



h 2 



X + —W 

go 



'0 



K go 



(Pl))+0(e"), (26) 



where the angular brackets denote averages over the parameter distributions, and where 
the integral quantities X, Y, Z, W are defined by equations (|T2l) . 



Proof: In the limit £k,Pk 0, the matrix element h — > ho- Since \ho\ < 1, the solutions are 
stable [22]. For nonzero fluctuations, the transfer matrix ((3]) takes the form of an elliptical 
rotation matrix [4], which can be written in the form 



M k = S(9 k ,L k ) 



cos 9k 
[1/L k ) sin 9 k 



-L k sin 9 k 

cos 9 k 



(27) 



where cos9 k = h k and L k = (sin 9 k )/g k . Theorem 4 of Ref. [4] shows that the growth rate 
for multiplication of elliptical rotation matrices of the form (12?]) . and hence the growth 
rate for the CC Random Hill's equation under consideration, can be written 



7 = -(sin 2 9 k )(r)l) + 0(e 



(2f 



where the perturbation r\ k is defined by L k = L (l + rj k ). 
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The length parameter L k of the elliptical rotation is given by 

Lk 



(29) 



9k 

Using the expressions for the matrix elements from equation (fTTjl in the definition (1291) of 
L k , and keeping only the leading order terms in e, we find 



Lk 



:i-^) 1/2 

#0 



1 + 



/,n [£ k X + p k Y] + -[£ k W + p k Z] 
'o 9o 



hi 



O e 



(30) 



where X,Y,Z,W are the integral quantities defined by equations (fl2l) or (I25I) . We can 
thus write the length parameter L k in the form L k = Lq(1 + r] k ), where 



(i - hiyi* 



go 



and 



[4x + + - [4W + . (3i) 



h 



o 



,9o 



Since we have defined the parameters l k and p k (which vary from cycle to cycle) to have 
zero mean, it follows that (rj k ) = 0. The quantity that appears in the growth rate from 
equation ( 128]) is the expectation value (sin 2 6 k ) (r) k ). The perturbation r\ k = 0(e), whereas 
to leading order sin 2 6 k = 1 — h\ + 0(e). As a result, the growth rate can be written in the 
form 

1 \ 2v 



7 



hi) 



T2 F + - Z 
«0 ^0 



+ O (e 3 ) . (32) 



The random variables £ k and are independently distributed, so that averages over the 
distributions can be rewritten in the form 



ho 



I -hi 



X + —W 

9o 



+ Pk 



h 



1 - 
h 



h 2 
a 



Y + —Z 

9o 



h 2 
a 



X + —W 



ho 



-l 2 



h 2 
a 



Y + —Z 

9o 



(Pt) ■ (33) 



I - nr, go 

The growth rate 7 thus has the form given by equation 

Theorem 2.3: Consider a Stochastic Hill's equation of the form ([6]) and a CC Random 
Hill's equation of the form (j2J), where the parameters have the forms A + £ k and q + p k . 
Both the parameter variations (£ k ,p k ) and the stochastic process £ are 0(e). To leading 
order in e, the growth rates of the solutions of the two types of differential equations are 
equivalent when the parameter variations (£ k ,p k ) are chosen according to 



and 



Pk 



(34) 
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where the integrals Ij and Jj are defined by equation f|24l) . and where the Sj are integrals 
of the stochastic process defined by 

Si = / 2/oi£^ and S 2 = / 2/02^ . (35) 
jo Jo 

The growth rates 7 are given by equation ( |26l ). 

Proof: The growth rates of both types of differential equations can be written in terms of 
the growth rates of the transfer matrix (under matrix multiplication). The growth rates 
will be the same if the matrix elements are the same (see equation Q). 

The Stochastic Hill's equation §6§ has a solution of the form given by equation ([9]), 
where yo(t) is the solution to equation ([6]) for the case where the right hand side vanishes. 
We consider the stochastic process £ and the correction y k to be O(e). To consistent order, 
y k (t) must obey the equation 

d 2 y k 



dt 2 



+ [\ + qQ(t)]y k = £, (36) 



If we multiply the above equation by yo{t) and then integrate over one period, from t = 
to t = Ti, we obtain the result 

/*7T 

yo(^)i/M -yMvki^) = / yoZdt. (37) 

Jo 

Now compare this result to that obtained from the CC Random Hill's equation where the 
parameters vary from cycle to cycle, 

y (n)y k (n) - y (n)y k (n) + 4 £ y 2 dt + p k jT y 2 Qdt = O (e 2 ) . (38) 

The matrix elements of the transfer matrices will be the same if equations (1371) and (J38l) 
are equal. But we obtain equations of these forms for both principle solutions, so this 
requirement results in two coupled equations 

+ PkJi = Si and £ k I 2 + Vk-h = S 2 , (39) 

where where Ij and Jj are the moments defined in equation f )24|) and the Sj are moments 
of the stochastic process defined through equations (|35|) . 

Equations (1391 provide two equations for two unknowns, and can be solved to find the 
perturbative quantities £ k and p k . We thus find the constraint of equation as claimed. 
The moment integrals Ij and Jj are defined by the solutions to the classical Hill's equation 
and are the same for all cycles. On the other hand, the moments Sj of the stochastic 
process will, in general, have different values for each realization (over each cycle labeled 
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by the index k). For a given cycle, equations fl34|) thus define the values of the 4 and p k 
that make the matrix elements for the CC Random Hill's equation that same as those for 
the Stochastic Hill's equation. This procedure thus results in distributions for the 4 and 
Pk- Finally with the distributions of the (equivalent) 4 and pk specified, the growth rate 
7 is given by equation (|26|) . • 

Corollary: Consider a stochastic Hill's equation in the alternate form 

% + [A + (q + £) Q(t)]y = . (40) 

The solutions to this differential equation will have the same growth rate as that of the 
CC Random Hill's equation when the parameter variations (£k,Pk) are chosen according 
to the constraints of equation (|34"|) . where the moments Sj have the alternate form 

Si = - [ W QyUdt and ~ 2 = - r QyUdt . (41) 
Jo Jo 

Proof: The proof is analogous to that of the previous result. • 

Remark: In the limit where the stochastic process £ has a correlation time r c — > 0, 
the integrals Sj in the numerators of equations (|34"|) vanish. As a result, the equivalent 
perturbations vanish, 

lim 4 = and lim p k = . (42) 

In this limit, the growth rate of either differential equation will be the same as that of the 
corresponding Classical Hill's equation. 

III. HILL'S EQUATION WITH SMALL FORCING PARAMETERS 

This section considers the CC Random Hill's equation fl2]) in the limit where the natural 
oscillation frequency A is constant and the forcing parameter g& = 0(e) where e < 1. This 
application is thus a special case of the general problem considered in the previous section. 

We first find solutions for the matrix elements h = yi(n) and g = yi(7r). The solutions 
can be expanded in orders of q k = 0(e), where the zeroeth order solutions are given by 

yoi(t) = cosv^At, and y 2 (t) = —j= sin \f\t . (43) 

v A 

The first order parts of the solutions obey the equation 

+ X V" + 9*^o = 0(e 2 ) , (44) 

where the parameter q k = 0(e). Next we multiply by yo and integrate over the interval 
[0,7r]. After integrating by parts twice, the expression becomes 

+ q k r Qyldt = 0(e 2 ) . (45) 



yoyk - yoyk 
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Note that this form holds for the perturbations for both the first and second principal 
solutions. Since the zeroeth order solutions (given by equation (143]) ) satisfy the boundary 
conditions, the perturbations y k and their derivatives y k must vanish at t — 0. The solution 
thus becomes 

yoj(n)y kj (ii) - y j(rc)y kj (Tc) + q k Jj = , (46) 

where the index j = 1,2 determines the principal solution and the moment integrals Jj are 
defined by equation ( |24"1) . The matrix elements h and g are given by 



h = y 1 (ir) = cos(p + y kl (ir) and g = y^n) = -VJ sirup + y kl (n) , (47) 

where we have defined tp = y/Xir to simplify the zeroeth order parts. Using these forms in 
the two equations (|46|) . we find two expressions for the two unknown matrix elements h 
and g: 

g cos <p + hVXsimp + q^ = C(e 2 ), (48) 

and 

h^l - L cos ip + q k J 2 = 0(e 2 ) . (49) 

VA g 

The matrix elements are thus given by 

h = cos ip -j= sinyj + 0(e 2 ) , (50) 

2y A 



and 

9 



-VXsmip + ^^f-sirSip-J^+Oie 2 ). (51) 
cosip \2 J 

Note that the expression for h is the same as the full solution for h in the limiting case where 
the barrier shape Q is a delta function [3]. In this same limit, the integral J\ = cos 2 ((^/2), 
and the expression for g also gives the exact result. Notice also that the results of Theorem 
2.1 reduce to equations (150!) and (15T|) in the limit £ k — > and q + p k — >■ q k (where equations 
[45] specify the zeroeth order solutions). 

Note that when cos<£> — > 0, the expression for g must be evaluated carefully. Both the 
numerator and the denominator vanish, and the limit must be evaluated. These values 
occur when A = (k + 1/2) 2 , where k G Z. 

Figure [TJ illustrates the fidelity of the first order approximation scheme. The figure 
shows the matrix element h as a function of A using the exact (numerically determined) 
form and the first order expression from equation (|50"1) . For this illustration, the parameter 
q k has constant value q k = 1/2. Note that the first order expression provides a good 
estimate. Smaller q k values will result in an even better approximation. 
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Figure 1: Matrix element h as a function of natural oscillation frequency A for Hill's 
equation with small forcing parameter = 1/2. Solid curve shows numerically determined 
result. Dotted curve shows the result for the approximation of equation ( 150|) . which is 
correct to leading order in q^. Dashed curve shows the result for the approximation of 
equation ( 1521 . which includes some higher order terms (see text). 
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For completeness, we note the following: Although the approximation (from equation 
15]) ) ignores terms that are second order in q k , after we make that approximation we 
can find the corresponding solution keeping all terms, i.e., we don't have to get rid of the 
remaining second order terms. We thus find the generalized solutions 



h 



Qk 



and 



9 



sin (p ± 



2 Qk 2 2 t t 

cos f + —sm f - q k J x J 2 



1/2 



cosy? 



—QkJi + sin 2 cp =F -\/A sin if 



2 Qk 2 2 t t 

cos + — sm f - q k JiJ 2 



1/2- 



(52) 



(53) 



where the signs are chosen so that the expressions reduce to those of equations ( 1501) and 
( IBTT) when the q\ terms are ignored. 

We now consider the growth rates. The transformation matrix of equation ([3]) can 
be written in the form of an elliptical rotation matrix, as in equation (12?]) . The length 
parameter of the rotation can be written in the form L k = Lq(1 + i] k ), where the T] k vary 
from cycle to cycle. For the case of symmetric variations, where (r] k ) = 0, the growth rate 
takes the form 



7 = !og 



(54) 



where 9k is the angle of the elliptical rotation matrix (see [4] and Theorem 2.2). 

Using the matrix elements found above, and dropping the subscripts on the h and g, 
we find 



cos 9 k = h = cos f — 



Qk 
2V\ 



smf + 0(e 2 



so that 



and hence 



9 k = f + 



qk 
2V\ 



sin 9i 



+ 0{t 



;i _ /^i/a = sinip + J* cosip + (e 2 ) 

2v A 



(55) 
(56) 
(57) 



The length parameter L k of the elliptical rotation matrix can then be written in the form 



sin 9i 



i + 



1 - 2Ji q k 
cos f sin f 2a/A 



1 



[l + Vk } + 0(e 2 



(5f 



where the final equality in equation (|58|) defines the perturbation r] k . Combining with the 
above results we find the following Theorem: 
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Theorem 3.1: The growth rate for the CC Random Hill's equation in the limit of constant 
A and small q = qu = 0(e) is given by 



7 = log 



]\ V cos (p 



(qt) (2J x -lV , ^ 2 



iX V cos ip 



+ 0(e 2 ) . (59) 



Note that this form is valid for the case of symmetric variations of the length parameter, 
i.e., where the distribution of the rjk is symmetric w.r.t. zero. 

Notice also that the expression for the growth rate is not defined when cos<^ — > 0. We 
can define the quantity J such that 

2Ji — 1 1 r 77 i— - 

J = — = / cos(2VXt)Qdt . (60) 



cos ip cos(VXn) 

When cosy? — > 0, \f~X — > (2k + l)/2, and both the numerator and denominator of the above 
expression vanish. To evaluate the quantity J, we can use L'Hopital's rule, which implies 
that 

J = ±— ft sin[(2A; + l)t\Qdt , (61) 

7T JO 

where the +(— ) sign arises for k even(odd). 

Figure H] illustrates how well this approximation scheme works. For the sake of definite- 
ness, the q^ are allowed to vary over a uniform distribution with amplitude A q . We expect 
the above results for the growth rate to be exact in the limit of small fluctuations, i.e., 
where the fluctuation amplitude A q — > 0. The figure shows the growth rates as a function of 
the amplitude A q for constant A = 1/2. Here the matrix elements are calculated using the 
first order expressions from equations (|50|) and (15 ip . The solid curve shows an estimate of 
the growth rate determined directly from matrix multiplication; note that the growth rate 
converges very slowly at small amplitudes, so the plot contains errors due to incomplete 
sampling. The dotted curve shows the growth rate calculated from equation (I54"|) . where 
the variations (77^) are found by numerical sampling. The dashed curve shows the growth 
rate from equation (|59p . which represents the main result of this section. Finally, even 
though the expressions derived here are correct only to first order in qj,, we can include 
the second order terms (in qk) for purposes of determining the rjk and hence the growth 
rate; for this case, the result is shown as the dot-dashed curve. All four of the curves are 
coincident for sufficiently small amplitudes A q . The analytic expression of equation (I5^|) 
represents the crudest approximation (dashed curve), and converges the slowest. 

Corollary: Consider a CC Random Hill's equation in the limit of small (but finite) forcing 
strength \qk\ -C 1, with constant oscillation parameter A, and where Q{t) is a function. In 
the limit A — > 00, the growth rate 7 — > 0. 
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Proof: In this limit, the growth rate is given by equation f l59|) . For nonzero fluctuations 
{qD > 0; the growth rate vanishes if and only if the integral Ji = 1/2, i.e., 

J x = dtQ(t) cos 2 (V\t) = 1/2. (62) 
Jo 

Using trigonometric identities, this expression can be written in the form 

I" dtQit) cos(2v / At) = 0. (63) 
J o 

In the limit A — > oo, this integral vanishes, and hence 7 — > 0. • 

Remark: The above argument works provided that Q(t) is a function. For the limiting 
case where Q is a delta-function (and hence a distribution), the integral of equation (1631 
does not necessarily vanish; nonetheless, the growth rate 7 oc 1/A in this limit [3], so that 
7 —¥ as A — > 00. 

For example, let Q = (2/n) sm 2 t. Then the quantity 2Ji — 1 becomes 

2Ji-l = t^^, (64) 

2vr v / A(A - 1) 

which vanishes as A — > 00. 

The general trend of decreasing growth rate 7 with increasing A is shown in Figure [3] 
for three choices for the barrier function. The various curves show the growth rates for 
Q(t) = (8/37r)sin 4 t (solid), Q(t) = (2/7r)sin 2 t (dashed), and Q(t) = 5([t] -tt/2) (dotted), 
where we have included the normalization constants, and where [t] denotes that the variable 
t is to be evaluated mod-7r. The amplitude of the q% fluctuations are the same for all three 
cases. 

IV. CONCLUSION 

This paper has generalized the Classical Hill's equation (0Q) to include random elements, 
where this treatment focuses on the case of small fluctuations with amplitude 0(e). This 
generalization can take two distinct, but related, forms: The parameters of the CC Random 
Hill's equation (j2J) vary from cycle to cycle, whereas the Stochastic Hill's equation (jSJ) or 
(JZJ) includes a stochastic process. 

Theorem 2.3 shows that the growth rates of these two types of (generalized) Hill's 
equations are the same when (for a given stochastic process £) the parameter variations 
of the CC Random Hill's equation are chosen according to equation ( 134"|) . Theorem 2.2 
provides an approximation for the growth rate, correct to 0(e 2 ), for the CC Random Hill's 
equation (and hence for the equivalent Stochastic Hill's equation). These growth rates 
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Figure 2: Growth rates for random Hill's equation in the limit of small random forcing 
parameter qk, plotted as a function of the amplitude A q of the fluctuations (for A = 1/2). 
The solid curve shows the result from direct matrix multiplication; dotted curve shows the 
result from equation (1541) . where the rjk are sampled numerically; the dashed curve shows 
the result from equation ( 1591) : the dot-dashed curve shows the growth rate calculated by 
including higher order terms (in qk). 
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Figure 3: Growth rates for random Hill's equation as a function of A, where the fluctuation 
amplitude (q%) = 1. The three curves show the result for Q = S([t] — n/2) (dotted), 
Q = (3/8?r) sin 4 t (solid), and for Q = (2/tt) sin 2 t (dashed). 
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7 ~ (pI) are nonzero except for isolated cases (see Figure [3]). As a result, for the CC 
Random Hill's equation, the (generalized) plane of parameters has no regions of stability. 
This result is not due to the fluctuations moving the matrix elements into the unstable 
regions of parameter space. One obtains nonzero growth rates even when all of the matrix 
elements are classically stable. In contrast, the Stochastic Hill's equation can be stable, 
where stability requires that the moments Ej of the stochastic process vanish (see equations 
[35] and [H]). Stability occurs in the limit where the correlation time r c — > 0. 

The results of this paper can be used in a wide variety of applications. As one example, 
orbit instabilities in extended mass distributions, such as dark matter halos, lead to a 
CC Random Hill's equation, as shown in Appendix A (see also Refs. [5,12]). As another 
example, the reheating problem at the end of inflation naturally leads to a Stochastic Hill's 
equation, as shown in Appendix B (see also Refs. [17,18,27]), and the growth rates can be 
calculated using Theorems 2.2 and 2.3. 
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APPENDIX A: 

RANDOM HILL'S EQUATION FROM ASTROPHYSICAL ORBITS 

Orbits in extended mass distributions are subject to an instability [5] that can be 
described by a CC Random Hill's equation (j2j). As one example, the density profile p{vj) 
of a dark matter halo has the form 

p{w) = po ^ , (Al) 
w 

where po is a density scale and the variable w is written in terms of Cartesian (x, y, z) 
coordinates through the relation 

^2 x 2 y 2 z 2 

w =^ + V 2+ V 2 ' (A2) 

where (without loss of generality) a > b > c > 0. The density field is thus constant on 
ellipsoids. The function F{w) approaches unity in the limit w — > so that the density 
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profile approaches the form p ~ 1/w. For this regime, one can find analytic forms for both 
the potential and the force terms [5]. 

When an orbiting body is initially confined to any of the three principal planes, the 
motion can be unstable to perturbations in the perpendicular direction. Consider an orbit 
initially confined to the x — z plane, with a small perturbation in the perpendicular y 
direction. In the limit \y\ <C 1, the equation of motion for the ^-coordinate takes the form 

^l + nlv = where & = - Jtt . (A3) 

dt 2 yy y ^ c 2 x 2 + a 2 z 2 + bVx 2 + z 2 V ' 

In this context, the time evolution of the coordinates (x, z) is determined by the original 
orbit and the motion is nearly periodic. As a result, the [x(t),z(t)} dependence of the 
parameter Q 2 provides a (nearly) periodic forcing term. The orbit has outer turning points 
which define a minimum value for Q 2 , which in turn defines the natural oscillation frequency 
Afc. The function Q 2 defined above can be written in the form 



y 
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n 2 = X k + Q k (t) } (A4) 

where the index k counts the number of orbit crossings, and the chaotic orbit in the original 
plane leads to different Xk and Quit) for each crossing. The shape of the functions Qk{t) 
are nearly the same, so that one can write Qk{t) = QkQ{t), where the forcing strength 
parameters g& vary from cycle to cycle. These forcing strengths qk are determined by the 
inner turning points of the orbit (weighted by the axis parameters [a, b, c]). Given the 
expansion of equation (IA4I) . the equation of motion ( 1A3I) for the perpendicular coordinate 
becomes a CC Random Hill's equation, with the form of equation (J2]). 



APPENDIX B: 

STOCHASTIC HILL'S EQUATION FROM REHEATING IN INFLATION 

In the inflationary universe paradigm [14] , the accelerated expansion of the universe is 
driven by the vacuum energy associated with a scalar field ip (or fields). During the phase 
of accelerated expansion, the energy density of the universe itself decreases exponentially 
and the cosmos becomes increasingly empty. This epoch is thought to take place when the 
universe is extremely young, with typical time scales of ~ 10~ 36 sec [19]. In order for the 
inflationary epoch to solve the cosmological issues it was designed to alleviate, the end of 
inflation must include a mechanism to refill the universe with energy [19]. This process is 
called reheating or preheating. 

During the reheating epoch, the equation of motion for the inflation field displays 
oscillatory behavior about the minimum of its potential. In order for the universe to 
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become filled with energy (reheat), the inflation field (p must couple to matter or radiation 
fields. One simple type of interaction arises from a coupling term in the Lagrangian of the 
form 



where x is another scalar field that represents matter (or radiation) and where the coupling 
constant g sets the interaction strength. The field x can be expanded in terms of its Fourier 
modes xt since these quantities evolve independently. The resulting equation of motion for 
the matter field modes xe then takes the form 



where Qe(t) is a periodic function (given by oscillatory behavior of the inflation field) and £ 
is a noise term that can be described by a stochastic process. Keep in mind that the index 
I refers to the Fourier mode. In the absence of fluctuations, the modes xe of the matter 
fields obey a Classical Hill's equation (JT]), which can be subject to parametric instability 
[17,18]. However, the noise perturbations £ [17,18,27] convert the reheating equation C IB2[) 
into a stochastic Hill's equation. As shown by Theorem 2.3, the growth rates for this 
stochastic Hill's equation are the same as an equivalent CC Random Hill's equation, where 
the conditions for equivalence are given by equation ( 155|) ; the growth rates can thus be 
calculated according to Theorem 2.2. 



£ in t = gpx 2 , 



(Bi) 




(B2) 



dt 2 
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